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REEFER  calculation  of  Jb(r)  and  0(r)  at  the  anode 
for  2.5  mm  cathode  dish  depth  and  10.0  mm  A-K  gap 
at  Vd  * 1.14  MV.  As  in  Figure  32,  Jb  has  been 
corrected  for  bipolar  flow,  whereas  $ is  from  the 
monopolar  calculation. 
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to  comparison  with  the  measured  profiles,  because  of  the 
current  collectors'  half-cm  diameter.  The  effect  is  evi- 
dent in  Figures  11(e)  and  (f ) , but  is  absent  in  (b) ; the 
profiles  in  parts  (a) , (c)  and  (d)  of  the  figure  do  not 
contain  measurements  near  the  edge  of  the  beam.  In  fact 
it  is  possible  that  in  part  (b) , the  enhancement  could 
have  fallen  between  the  two  outermost  radii  of  measurement, 
which  were  slightly  over  1 cm  apart. 

• REEFER  predicts  a central  peak  in  the  J,  profile,  as 
observed,  but  the  calculated  peaks  are  smaller  than  in 
the  measured  profiles.  The  peak  is  most  probably  caused 
by  beam  pinch  and/or  beam  compression,  and  must  occur  if 
electron  trajectory  angles  in  the  r-z  plane  go  zero 
slower  than  linearly  with  decreasing  r.  Reasons  for  the 
calculated  peaks  being  smaller  than  observed  are  (1) 

REEFER  could  not  accommodate  a converging  guide  field, 
so  a paraxial  field  had  to  be  used  in  the  calculations 
(which  will  not  be  necessary  in  the  future,  as  the  code 
is  to  be  modified  to  include  the  radial  magnetic  field 
component) ;and  (2)  REEFER  computed  the  electron  trajectories 
in  the  monopolar-mode  diode,  so  the  amount  of  pinching  is 
less  than  would  occur  in  the  higher-current  bipolar  mode. 

• The  two  points  just  mentioned  may  also  explain  why  the 
calculated  beam  radius  at  the  anode  (4.9  cm  in  Figure  30) 
is  non-negligibly  larger  than  the  measured  beam  radius 
(typically  4.2  cm  in  Figure  11).  Furthermore,  the  J.(r) 
profile  was  measured  1 cm  from  the  anode,  where  guide 
field  flux  tubes  were  5 percent  smaller  in  radius  relative 
to  the  anode  plane.  We  note  also  that  the  damage  pattern 
radius  could  be  up  to  two  gyroradii  (or  1 cm)  larger  than 
would  be  the  case  for  a beam  with  zero  transverse  energy. 
Because  of  all  these  effects,  direct  comparison  of  the 
calculated  and  measured  beam  radii  (Figures  11  and  30) 
does  not  have  an  immediately  obvious  interpretation. 


• Azimuthal  asymmetries  in  the  measured  J.  profiles  in 
Figure  11  make  them  difficult  to  compare  against  the 
calculated  profiles,  as  regards  details  of  shape  between 
the  central  peak  and  the  edge  enhancement.  The  measured 
profiles  seem  generally  more  "hilly,"  though  of  the  same 
type  of  shape  as  the  calculated  profiles.  Quantitatively, 
J,  amplitudes  are  in  reasonable  agreement  in  most  cases; 
appreciable  discrepancy  occurs  in  the  interior  of  the 
beam  with  the  deeper-dished  cathode,  especially  in  Fig- 
ure 11(e)  at  radii  of  about  1 and  2%  cm.  Based  on  the 
measured  J.  values  at  those  locations,  anode  dose  would 
have  been  marginally  below  the  bipolar  flow  threshold 
(locally)  at  the  time  shown;  the  monopolar  REEFER  values 
of  J.  do  agree  with  the  measurements.  At  larger  radii 
(where  the  A-K  was  smaller  by  2 to  2*s  mm)  , anode  dose 
should  have  been  high  enough  to  trigger  bipolar  flow  by 
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the  time  shown.  This  case  suggests  that  the  cathode 
contour  can  be  used  to  localize  the  region  of  bipolar 
flow.  Referring  again  to  Figure  11(e),  the  naturally- 
occurring  central  peak  in  the  monopolar  profile  must 
have  been  high  enough  to  trigger  bipolar  Dflow  at  the  axis. 

_2 

• The  theoretical  “ d dependence  of  paraxial  space- 
charge-limited  flow  between  planar  electrodes  is  evident 
in  the  REEFER  results.  The  REEFER  A-K  gap  at  the  axis 
was  11.15  mm  for  the  1.5  mm-deep  cathode,  and  12.5  mm 
for  the  2.5  mm-deep  cathode,  giving  a d2  ratio  of 

(11 . 15/12 . 5 ) 2 = 0.796.  The  calculated  values  in  the 
vicinity  of  the  axis  (averaged  over  the  region  r < 1.5  cm) 
are  3.2  kA/cm2  and  3.9  kA/cm2  for  the  2.5  mm  dish  and  1.5  mm 
dish,  respectively;  3. 2/3. 9 = 0.82,  in  good  agreement  with 
the  d'2  ratio. 

• As  calculated  by  REEFER,  the  electron  angles  at  the  anode 
plane  are  generally  smaller  than  predicted  by  the  net 
magnetic  field-following  model  (Section  4).  As  an  example, 
we  consider  the  REEFER  <J>(r)  results  shown  in  Figure  32;  at 
radii  of  2.35  cm,  3.56  cm,  and  4.9  cm, the  field-following 
model*  gives  <J> f f / 4>  = 1.28,  1.30,  and  1.66,  respectively. 

The  field-following  model  neglects  the  diode's  axial  elec- 
tric field,  which  brings  electron  trajectories  straighter 
toward  the  anode  than  would  pure  magnetic  f ield*- following 
motion.  Furthermore,  Figure  30  shows  non-negligible  radial 
components  of  electron  velocity  at  the  anode  due  to  beam 
self-pinch.  The  effect  of  the  radial  velocity  component 

is  to  increase  <p ; this  effect  (included  in  the  4>  ( r ) plot 
in  Figures  32  and  33)  also  is  neglected  in  the  field- 
following model. 

6.1.3  Limitations  of  REEFER.  Several  inadequacies  were 
encountered  in  REEFER  in  the  course  of  the  studies  just  reported. 
Three  points  have  already  been  mentioned:  (1)  numerical  insta- 
bility of  regions  with  strongly  non-uniform  current  density,  such 
as  pinched  flow  regions,  due  to  coarse  field-updating  procedure; 

(2)  inflexibility  of  electrode  contouring  due  to  a Ar  = Az  re- 
quirement in  the  code's  Poisson-solving  algorithm;  (3)  REEFER 
can  not  simulate  bipolar  flow  or  re-entrant  electrons,  which 

*with  B = 0,  as  in  the  REEFER  calculation.  Also,  both  4>,f  and 
the  REEFER  4>-values  are  from  the  monopolar  flow  case. 
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we  have  shown  to  be  important  features  of  diode  behavior  in  the 
experimental  domain  of  interest.  To  these  points  we  add  that 
REEFER  is  an  equilibrium  code,  incapable  of  accounting  for  the 
effects  of  time-varying  diode  voltage  on  particle  motion.  This 
may  not  be  a serious  limitation  for  monopolar  (electron)  diodes 
but  could  be  important  for  bipolar  diode  simulation,  because  of 
relatively  long  ion  transit  times. 


6.2  BEAM  TRANSPORT  AND  COMPRESSION  STUDIES  USING  THE  COMPUTER 
CODE  REFLEX 


The  computer  code  REFLEX  was  developed  in  this  program  to 
numerically  evaluate  the  Young-Spence  model  of  beam  transport  in 
a converging  Bz  guide  field  (Reference  3,  reproduced  in  the 
appendix  of  this  report) . The  model  was  originally  motivated  by 
current  transport  efficiency  measurements  (see  Reference  3) 
which  were  higher  than  could  be  explained  by  calculations  of 
simple  adiabatic  motion  of  beam  electrons.  The  adiabatic  calcu- 
lations predicted  too  much  mirroring  (reflection)  of  electrons 
in  the  spatially  increasing  field.  The  Young-Spence  model  in- 
voked electron  axial  reflexing  between  the  Bz  mirror  (or  lens) 
field  and  the  cathode.  Reflexing  causes  untransmitted  electrons 
to  repeatedly  re-encounter  the  anode,  where  scattering  produces 

a non-adiabatic  change  in  pitch  angle.  If  an  electron's  new 
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angle  is  below  the  "loss-cone"  angle  sin  (M  ) , where  M is  the 
ratio  of  peak  Bz  to  Bz  at  the  anode,  the  electron  will  then 
propagate  through  the  lens.  In  our  case,  the  loss-cone  angle 
was  33  degrees.  Since  a radial  "walk"  accompanies  scattering, 
an  electron  can  leave  the  system  (radially)  before  scattering 
puts  its  angle  below  the  critical  value  for  transmission.  The 
Young-Spence  model  computes  transport  efficiency  through  the 
lens  as  the  outcome  of  these  two  competing  effects. 


A detailed  numerical  treatment  of  the  model  was  desirable 
because  its  analytic  version  (Reference  3)  cannot  take  account 
of  the  variation  of  injected  electron  energy  with  time  or  the 
variation  of  injected  beam  current  density  with  radius.  Further- 
more, the  analytic  model  computes  only  charg-e  transport  efficiency, 
and  unlike  the  new  REFLEX  code,  cannot  calculate  other  important 
features  of  the  transmitted  beam,  such  as  total  current  versus 
time,  current  density  profile  versus  time,  and  radius-  and  time- 
resolved  energy  and  angle  distributions.  These  transmitted  beam 
quantities  Will  be  important  as  input  data  for  Monte  Carlo  code 
computations  of  the  energy  deposition  profile.  Present  emphasis 
is  on  refining  the  physics  of  the  code  based  on  experimental 
feedback. 

As  mentioned  in  the  Introduction  (Section  1) , it  initially 
seemed  desirable  to  couple  the  transport  code  to  REEFER,  i.e., 
to  use  REEFER  diode  calculations  as  input  to  the  transport  calcu- 
lation. The  software  for  interfacing  the  two  codes  was  in  fact 
developed,  but  never  used,  because  the  injected  electron  angles 
computed  by  REEFER  were  for  the  monopolar  flow  mode  and  so  were 
unrealistically  low  as  it  turned  out.  With  the  converging  guide- 
field  geometry,  realistic  electron  angles  of  injection  were  felt 
to  be  essential  for  a meaningful  test  of  the  transport  code.  We 
therefore  employed  the  field-following  model  as  the  best  way  of 
estimating  angles  of  injection;  this  simple  model  was  found  to 
be  consistent  with  angles  measured  near  the  anode  (Section  4) . 

The  injected  beam  was  also  taken  to  have  uniform  current  density, 
a reasonable  approximation  of  the  data  l Figure  11(a)  and  ( b ) ] for 
the  purpose  of  initial  tests  of  the  code. 
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The  new  code,  REFLEX,  is  phenomenological  in  nature.  It 
segments  the  injected  beam  in  radius  and  time  and  then  tracks 
the  evolving  transport  status  of  each  beam  segment  according 
to  analytic  laws  governing  electron  motion  and  scattering  (diffu- 
sion) . The  output  beam  is  fully  described  as  a time-varying 
electron  distribution  over  radius,  energy,  and  angle.  Certain 
adjustable  parameters  must  be  specified:  (a)  the  radius  at 
the  anode,  outside  of  which  electrons  are  assumed  to  be  lost 
from  the  system;  (b)  the  scattering  angle  probability  function 
describing  the  effect  of  each  encounter  with  the  anode  foil  and 
the  intra-diode  fields.  The  complex  structure  of  the  latter 
has  prevented  analytic  calculation  of  the  scattering  distribution, 
so  the  code  assumes  re-injection  with  angles  uniformly  distributed 
up  to  a maximum  angle  which  depends  on  electron  energy  and  incident 
angle.  Item  (a) , the  containment  radius  at  the  anode  can  be  identi- 
fied experimentally  in  two  ways.  For  the  present  study,  it  was 
determined  by  tracing  back  the  guide-field  flux  tube  at  the  edge 
of  the  compressed-beam  damage  pattern.  The  other  way  would  be 
to  place  various  plates  with  apertures  of  different  diameters 
against  the  downstream  side  of  the  anode,  in  turn  (with  each 
plate  thick  enough  to  absorb  incident  electrons) . The  contain- 
ment radius  would  be  defined  by  the  smallest  aperture  which  does 
not  affect  beam  transport. 

6.2.1  REFLEX  Calculations  Compared  with  Experimental  Results. 


The  waveforms  of  injected  electron  current  and  energy  that  were 
used  as  inputs  to  the  REFLEX  calculations  were  the  diode  waveforms 
of  Pulse  4048  (Figure  4b).  As  already  described,  the  injected 
radial  current  profile  was  assumed  uniform,  with  electron  angles 
given  by  the  field-following  model  in  the  diode.  The  target  plane 
for  which  beam  transport  was  calculated  was  taken  to  be  the  focal 
plane  of  3:1  Bz  lens,  with  the  same  field  strength  at  the  anode  as 
in  our  experiments;  the  anode  was  also  as  in  the  experiments 
(^-mil  titanium) . Because  the  guide  field  is  not  subject  to  radial 
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dilation  in  the  code  (as  evidently  it  was  experimentally)  the 
calculated  transport  efficiency  and  radial  current  profile  shape 
are  independent  of  target  position  at  all  points  downstream  from 
the  focal  plane.  Experimentally,  the  lens  ratio  was  3.3:1  before 
beam  injection,  with  the  target  (for  compressed  beams)  at  the 
3:1  position,  5 cm  downstream  from  the  lens  plane.  We  assume  in 
the  code  that  lens  dilation  during  the  pulse  reduced  the  field 
upstream  from  the  target,  so  that  the  actual  field  during  beam 
transport  increased  monotonically  from  anode  to  target,  where 
line-tying  maintained  the  3:1  value. 

Figures  34  through  37  present  REFLEX  calculations  of  the 
transported  and  compressed  beam.  The  figures  show  total  current 
transport,  total  energy  transport,  radial  current  profile,  and 
mean  electron  angle.  (Other  information  which  the  code  can  pro- 
vide, such  as  radius-  and  time-resolved  electron  angle  and  energy 
distributions  in  the  transported  beam,  were  not  available  in  time 
for  inclusion  in  this  report.) 

The  ratio  of  the  time-integrals  of  the  current  waveforms  in 
Figure  34  gives  a calculated  charge  transport  efficiency  of  80  per- 
cent, as  compared  with  the  experimental  value  of  95  percent  (see 
Section  2) . Total  energy  transport,  in  Figure  35,  is  calculated 
at  71  percent,  as  compared  with  an  estimated  86  percent  experi- 
mentally (based  on  the  measured  charge  transport,  as  described 
in  Section  2) . The  calculated  transport  values  are  about  83  per- 
cent of  the  experimental  values.  This  modest  but  non-negligible 
discrepancy  indicates  a need  for  further  refinement  of  the  code's 
modeling  of  non-adiabatic  scattering  when  reflexing  electrons  re- 
encounter the  anode  and  the  diode  fields.  It  is  encouraging  that 
the  transported  current  waveform  in  Figure  34  exhibits  the  same 
qualitative  behavior,  relative  to  the  injected  current,  as  does 
the  experimental  waveform  (Figure  5a  in  Section  2) : very  high 
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Figure  34  REFLEX  current  transport  calculation 
through  3:1  lens.  Injected  electron 
beam  current  (Ijnj)  and  transported 
current  (I)  versus  time.  Total  charae 
transport  efficiency  is  80  percent. 
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Figure  35  REFLEX  energy  transport  calculation 

through  3:1  lens.  Cumulative  injected 
beam  energy  (£^nj)  and  transported 
beam  energy  (£)  versus  time.  Final 
transport  efficiency  is  71  percent. 


current  transport  efficiency  for  the  first  20  ns,  followed  by 
85  ns  of  reduced  efficiency,  followed  by  higher  efficiency  for 
the  rest  of  the  pulse.  The  physical  reason  for  the  middle  period 
of  reduced  transport  efficiency  is  the  relatively  high  diode 
voltage  during  that  time.  Higher  electron  energy  increases  the 
electron  gyroradius,  which  in  turn  increases  radial  losses  asso- 
ciated with  anode  re-encounters  (Reference  3;  see  Appendix). 

Figure  36  shows  the  calculated  current  density  radial  pro- 
file at  the  time  of  peak  injected  current.  Recalling  that  the 
calculation  used  a uniform  (flat)  injected  current  profile,  the 
figure  represents  a theoretical  indication  that  beam  compression 
may  cause  considerable  profile  distortion.  Unfortunately,  the 
experimental  profile  measurements  at  the  compressed  beam  position 
(Section  3)  are  too  ambiguous  due  to  azimuthal  non-uniformities 
to  serve  as  a test  of  this  effect.  Physically,  the  calculated 
profile  distortion  is  a result  of  radial  diffusion  and  loss  of 
electrons,  brought  about  by  scattering  at  the  anode  (Reference  3). 
The  altered  profile  has  approximately  the  "correct"  amplitude  in 
its  central  region  (i.e.,  three  times  the  injected  current  density, 
corresponding  to  factor-of-three  area  compression) , but  is  de- 
ficient in  the  outer  portion  of  the  beam,  because  radial  escape 
occurs  from  the  beam  surface. 

The  average  electron  angle  plotted  versus  time  in  Figure  37 
is  the  simple  statistical  instantaneous  mean  value.  Different 
angle  distributions  can  have  the  same  mean  value,  so  the  latter 
is  not  the  same  as  the  instantaneous  effective  angle  for  electron 
energy  or  charge  deposition  in  a target  material.  Therefore, 
the  mean  angle  in  Figure  37  should  not  be  simply  compared  with 
the  measured  effective  angles  in  Section  4 (Figure  19).  Quali- 
tatively, the  time  dependence  of  the  average  angle  in  Figure  37 
differs  from  that  of  the  effective  angles  in  Figure  19,  in  that 
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Figure  36  REFLEX  calculation  of  compressed 
current  density  profile  (J)  at 
the  target  with  2.9:1  lens  ratio, 
at  the  time  of  peak  injected  current 
Maximum  J is  approximatel v 3.3  times 
the  injected  (uniform-profile)  value 
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Fiaure  37  REFLEX  calculation  of  average  anqle 
in  3:1  compressed  beam  versus  time. 
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the  latter  decrease  markedly  as  diode  voltage  falls  toward  the 
end  of  the  pulse,  whereas  the  calculated  average  angle  does  not. 
This  discrepancy  might  indicate  inaccuracies  in  the  field-follow- 
ing model  of  diode  electron  angles,  which  was  used  to  numerically 
generate  the  injected  electron  angles  for  REFLEX.  The  compressed- 
beam  electron  angle  calculation  would  best  be  tested  against  the 
data  by  using  the  full  calculated  angle  distribution  as  input  for 
a Monte  Carlo  code  computation  of  the  beam  current  transmission 
profile.  This  computed  profile  could  then  be  compared  directly 
with  filtered  Faraday  cup  data. 


SECTION  7 

SUMMARY  AND  RECOMMENDATIONS 

The  main  achievements,  observations  and  conclusions  of  the 
Electron  Beam  Diagnostics  Program  can  be  summarized  as  follows: 


1.  The  instantaneous  cross-sectional  profile  of  beam 
current  density  has  been  measured  in  350  kA,  1.2  MV, 

50  kJ  beams,  both  near  the  anode  and  after  factor-of- 
three  area  compression.  The  application  of  these  measure- 
ments to  high-dose  calorimetry  has  been  demonstrated. 

2.  Near  the  anode,  the  current  density  profile  data  are 
in  general  agreement  with  simulations  performed  using  the 
REEFER  diode  code,  provided  the  simulation  results  are 
adjusted  according  to  bipolar  flow  theory.  Experimental 
beam  non-uniformities  made  it  impossible  to  rigorously 
test  the  simulation. 

3.  The  quantitative  accuracy  of  filtered  Faraday  cup 
measurements  of  charge  deposition  profile  has  been 
established,  and  filtered  Faraday  cup  data  has  been 
analyzed  to  give  the  first  time-resolved  measurements 
of  electron  angles  in  intense  e-beams.  This  capability 
allows  the  accurate  measurement  of  high-dose  energy 
deposition  profiles  in  samples  of  practical  interest. 

4.  Beam  compression  by  a factor  of  3 in  area  was  ac- 
complished at  ~ 95  percent  charge  transport  efficiency 
and  ~ 85  percent  energy  transport  efficiency,  giving 
peak  dose  of  about  4000  cal/gm  in  carbon  targets. 

5.  A new  computer  code,  REFLEX,  has  been  developed  to 
numerically  simulate  the  Young-Spence  model  of  beam 
transport  and  compression  in  axial  magnetic  guide  fields. 
REFLEX  calculations  of  compressed-beam  transport  efficiency 
are  about  17  percent  less  than  measured  values,  suggesting 
a need  for  more  detailed  modeling  of  electron  scattering 

in  diode  re-encounters. 
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6.  Modification  of  diode  behavior  due  to  reflected 
electrons  has  been  experimentally  demonstrated,  docu- 
menting the  need  to  treat  this  effect  in  diode  simula- 
tion codes. 


Based  on  this  program's  findings,  and  in  view  of  the  overall 
program  goal  described  in  the  Introduction,  we  recommend  continu- 
ing efforts  in  the  following  areas: 


7 . 1 EXPERIMENTAL 


• Electrode  surface  treatment  tests,  aimed  at  subduing  diode 
impedance  collapse  mechanism (s)  in  order  to  achieve  beam 
azimuthal  uniformity,  detailed  control  of  the  radial  current 
profile  shape,  and  reproducibly  full-duration,  full-energy 
pulses . 

• Compressed  beam  transport  studies  (such  as  lens  ratio,  anode 
thickness,  and  cathode  contour  variations)  aimed  at  further 
revealing  the  mechanisms  controlling  beam  transport. 

• Radius-resolved  measurements  of  charge  deposition  profile 
and  electron  angle  by  the  filtered  Faraday  cup  technique. 

• Experimental  studies  to  define  the  domain  of  validity  of 
the  single-angle  beam  model. 


I 
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7 . 2 THEORETICAL/CALCULATIONAL 

• Diode  simulation  development,  including  the  effects  of 
bipolar  flow,  re-entrant  electrons,  spatially  varying 
guide  field,  and  time-varying  voltage. 

• Transport  theory  studies,  aimed  at  identifying  the  mechan- 
ism (s)  which  give  measured  transport  efficiency  in  excess 
of  presently  predicted  value. 

• Monte  Carlo  code  re-organization,  to  include  time-varying 
electron  angle. 

• Calculations  to  identify  single-angle  beam  model  range  of 
validity  to  complement/support  experimental  study. 

• Calculational  study  of  requirements  for  energy  deposition 
profile  radial  uniformity  (angle/f luence/dose/prof ile  shape 
trade-off) . 
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APPENDIX 

MODEL  OF  MAGNETIC  COMPRESSION  OF  RELATIVISTIC  ELECTRON  BEAMS 


Model  of  magnetic  compression  of  relativistic  electron 
beams 

T.  S.  T.  Young  and  P.  Spence 

Physics  International  Company.  San  Leandro.  California  945  77 
(Received  30  June  1976) 

A theoretical  model  for  the  magnetic  compression  of  relativistic  electron  beams  is  described  Calculations 
of  the  beam  fraction  transmitted  through  the  magnetic  mirror  are  in  agreement  with  several  experiments 
covenng  a broad  parameter  range 


PACS  numbers  84  80  Dz.  52.40.  Mj 

There  have  been  several  experimental  studies  on  the 
compression  of  relativistic  electron  beams  in  low- 
pressure  gas  using  externally  imposed  fields.  At  Sandia 
Laboratories,1  magnetic  compression  up  to  mirror  ra- 
tios in  excess  of  6 : 1 of  a low  v/y  beam  agreed  with  the 
adiabatic  prediction  that  the  beam  should  conserve  its 
magnetic  moment  while  propagating  along  the  field  lines. 
The  transmitted  beam  fraction  through  the  mirror  was 
simply  the  beam  population  in  the  magnetic  loss  cone. 

At  Cornell  University,  Davitian  and  Nation2  obtained 
rather  high  transmitted  beam  fractions  at  mirror  ratios 
of  3 :1.  In  addition  they  found  the  compression  of  the 
beam  cross-sectional  area  to  be  inversely  proportional 
to  the  mirror  ratio  Af,  indicating  that  the  beam  elec- 
trons were  tied  to  the  vacuum  field  lines.  In  experi- 
ments with  high  v/y  beams  at  Physics  International 
Company, 3-4  transmission  fractions  were  repeatedly 
observed  to  be  much  higher  than  the  adiabatic  predic- 
tion while  the  beam  cross-sectional  area  appeared  to 
scale  as  AC1.  Therefore,  a theory  is  needed  to  explain 
the  discrepancy  among  the  various  experiments  in  that 
the  adiabatic  theory  accounts  for  the  beam  transmission 
in  some  experiments  but  is  inadequate  for  the  others. 

It  is  proposed  here  that  the  beam  electrons  which  are 
adiabatically  reflected  at  the  magnetic  mirror  can  be 
trapped  between  the  mirror  and  the  electric  potential 
at  the  cathode.  The  trapped  electrons  nonadiabatically 
change  their  pitch  angles  as  well  as  their  guiding  flux 
tubes  when  passing  through  the  anode.  Hence,  some 
fraction  of  the  trapped  electrons  can  eventually  be 
transmitted  through  the  mirror  due  to  the  nonadiabatic 
change  of  their  pitch  angles,  while  the  remainder 
escape  radially  due  to  the  random  changes  of  their 
guiding  centers  Beam  transmission  is  enhanced  over 
the  simple  adiabatic  (single-oass)  prediction  when  the 
radial  escape  rate  is  slow  compared  with  the  pitch 
angle  scattering.  The  criterion  for  the  transmission 
enhancement  to  occur  is  found  to  be 

p*//Sr,(l-cos<’c),',<l,  (1) 

where  r,  is  the  beam  radius  at  the  anode  in  centimeters, 
8,  is  the  loss  cone  angle,  and  p,  (cm)  1.7(yJ-  l)1'- 
B,  (kG)  is  the  full  beam  gyroradius  in  the  plasma 
channel. 

We  assume  that  an  electron  beam  with  uniform  den- 
sity, which  is  generated  in  a vacuum  diode,  propagates 
in  an  axial  magnetic  guiding  field  and  is  compressed  in 
a mirror  field  at  the  end  of  the  transport  region.  The 
field  is  assumed  to  be  strong  enough  to  prevent  beam 
pinching  in  the  diode, 5 i.  e. , 
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B,  (kG)  7 0.U,  (kA)  y/d  (cm)  v.  (2) 

In  Eq.  (2),  l.  is  the  beam  current,  <1  is  the  cathode- 
anode  distance,  and  v/y  is  the  usual  Budker  parameter. 
In  the  transport  and  compression  region,  we  assume 
that  the  beam  is  fully  current  neutralized  and  that  the 
beam  electrons  are  adiabatically  guided  by  the  vacuum 
field  (B  - B,i  just  beyond  the  anode  and  I B I = MB,  ■ B, 
at  the  mirror).  Diamagnetic  effects  are  assumed  unim- 
portant because  in  most  of  the  experiments1"4  the  beam 
cross  sections  appear  to  compress  adiabatically.  The 
configuration  of  the  beam  environment  is  schematically 
plotted  in  Fig.  1 for  illustration. 

In  the  diode,  the  uniform  beam  produces  an  azimuthal 
magnetic  field,  B,  5,  which  adds  vectorially  to  the  ap- 
plied axial  B,.  For  the  cases  where  the  diode  electric 
field  is  small  compared  with  the  total  B field,  beam 
electrons  are  accelerated  by  the  electric  field  along 
the  total  B field  and,  at  the  same  lime,  gain  a small 
E x B drift  in  the  radial  direction.  The  electrons  then 
encounter  the  anode  foil  across  which  the  spiral  diode 
B field  changes  into  the  axial  B,.  Assuming  that  the 
transition  sheath  at  the  anode  foil  is  much  thinner  than 
a beam  gyroradius,  the  beam  electrons  streaming 
across  the  anode  sheath  into  the  plasma  channel  obtain 
a pitch  angle  e,  with  respect  to  B„  where 

8,  ^ (r/r,)  tan"1  (0. 2 /„  > , B,)  * 8,  r/r„  (3 ) 

Scattering  in  the  thin  anode  foil  induces  an  additional 
pitch- angle  spreading  centered  about  the  pitch  angle 
8r,  and  the  corresponding  energy  absorption  is  negligi- 
ble. The  angular  spreading  s,  (which  is  a function  of  r 


Vlt, . i.  schematic  configuration  of  the  diode,  beam  transport, 
and  beam  compression  regions. 
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because  the  beam  incident  angle  0r  and  hence  the  effec- 
tive anode  thickness  depend  on  r)  will  be  assumed  to  be 
constant  for  simplicity.  For  low  > beams  and  thin 
titanium  foil® 

«,*(20/y)K  (cm)|,'*1  (4) 

where  d,  is  the  anode  foil  thickness. 

We  assume  a two-dimensional  monoenergetic  beam 
distribution  function 


A»./>,r)  = 


erffir  «,) 


exp 


(5) 


where  />  is  the  electron  momentum,  6 is  the  usual  delta 
function,  and  erf  is  the  error  function.  The  correspond- 
ing beam  fraction  within  the  loss  cone  which  can  be 
adiabatically  transmitted  through  the  magnetic  mirror 
is 


«(»•)* X./  <10  St  dl>PA0,P,r) 

= {erf[(flc  + 0r)/$,]  + erf[(«c-  er)/e,\){2  erflx/d,)!'1, 


(6) 


where  «easin*1(t/Af,/2)  is  the  loss  cone  angle.  The 
total  adiabatic  transmitted  fraction  is  simply 


drrR(r) 


1 <2  ’V»)  / * 


x exp 


(7) 


The  analytic  expression  (7)  gives  almost  identical 
results  to  the  adiabatic  numerical  calculation  by- 
Freeman  and  Poukey. 1 

The  adiabatically  reflected  beam  electrons  then  turn 
back  toward  the  diode,  pass  through  the  anode  foil, 
bounce  back  from  the  cathode  potential  hill,  and  again 
pass  through  the  anode.  The  mean  squared  angular 
scattering  gained  by  the  beam  electrons  by  twice  cross- 
ing the  anode  foil  (thin  titanium  foil  assumed)  is 
roughly 

(SO1)  ai  (400/>J)rf,  (cm)  + 2«i,  (8) 

where  the  angular  brackets  denote  the  usual  ensemble 
average.  The  first  term  in  Eq.  (8)  is  due  to  the  thin 
foil  scattering.  The  second  term  is  due  to  the  nonadia- 
batic  crossings  of  the  anode  sheath.  One  may  substitute 
Eq.  (8)  for  the  angular  scattering,  Eq.  (4),  to  calculate 
the  beam  fraction  which,  after  reflection,  can  now  be 
adiabatically  transmitted.  However,  we  note  that  for 
many  experiments  the  value  of  <««*)  in  Eq.  (8)  is  of  the 
order  of  unity.  In  such  cases,  the  reflected  beam  par- 
ticles are  roughly  uniformly  distributed  in  pitch  angles. 
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Hence,  the  fraction  of  beam  scattered  into  the  loss 
cone  per  reflection  is  simply  l-cos0„  which  is  the  coni- 
cal surface  area  in  the  loss  cone  divided  by  the  total 
spherical  surface  area.  The  number  of  mirror  reflec- 
tions required  for  an  average  electron  to  be  scattered 
into  the  loss  cone  is 

N-  (1  - cos«e)*‘  a 2»;2.  (9) 

On  the  other  hand,  there  is  a radial  scattering  of  the 
electron  guiding  center  and  it  is  related  to  the  pitch- 
angle  scattering  by 

^ do) 

When  an  electron  is  scattered  to  a radius  comparable 
to  the  cathode  radius,  it  is  no  longer  trapped  by  the 
diode  electric  field.  The  number  of  mirror  reflections 
L required  for  an  electron  to  be  radially  lost  can  be 
roughly  defined  by  requiring  the  accumulative  (Sr2)  to 
be  equal  to  the  square  of  the  initial  particle  radial 
distance  from  the  beam  edge,  i.  e. , 

£ <8r*  > (11) 

<■1  • , 

where  r is  the  initial  radial  position  of  the  electron. 

When  (8S2)  is  of  the  order  of  unity  (i.  e. , 90°  scatter- 
ing), Eq.  (11)  can  be  rewritten 

L =*  (r,-  r)2/p£  a 2(r,-  r)\>p\  (12) 

by  using  Eq.  (10),  where  is  the  average  gyroradius 
of  the  electron.  An  electron  will  be  scattered  into  the 
loss  cone  and  transmitted  through  the  mirror  before  it 
is  radially  scattered  out  of  the  beam  channel  if  L * S, 
which  is  equivalent  to 

r/r,  *.  1 - p,/ v"5  r,(l  - cos8c)1/J  (IS) 

from  Eqs.  (9)  and  (12).  Equation  (13)  defines  a core 
area  at  the  anode  in  which  all  the  electrons  can  eventual- 
ly transmit  through  the  mirror,  while  outside  the  core 
area  only  the  adiabatic  fraction  [see  Eq.  (6) ) will  be 
transmitted.  Transmission  efficiency  is  enhanced  over 
the  adiabatic  prediction  in  Eq.  (7)  only  if  the  right-hand 
side  of  Eq.  (13)  is  positive. 


FIG.  2.  Comparison  between  calculated  transmission  (solid 
curve)  and  data  from  Ref.  4 for  various  mirror  ratios.  The 
triangles  and  circles  represent  20-  and  100-cm  transport  data, 
respectively.  The  corresponding  adiabatic  transmission  is  also 
plotted  (dashed  curve)  for  comparison. 
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TABLE  I.  Comparison  between  calculated  transmitted  fraction 
(T^)  with  three  experiments  (T.^,). 


References 

1 

2 - 

3 

/,  <kA) 

70 

40 

400 

y 

4.91 

1.78 

2.5 

v/y 

0.9 

1.4 

9.9 

rt  (cm) 

7.2 

3.0 

4.4 

B,  (kG) 

1.94 

2.67 

18.2 

B,  (kG) 

2.  5 

4.0 

9.0 

M 

3.7 

"TTN 

3 

4 

69  ±6 

41*5 

83-96 

75  * 5 

0.66 

/ 

0.59 

1.1 

0.28 

0.316 

0.52 

•c 

0.547 

0.40 

0.65 

0.524 

itt) 

65 

_ 

40 

83 

70 

We  now  proceed  to  compare  the  transmission  calcula- 
tion from  Eqs.  (6),  (7),  and  (13)  with  several  experi- 
ments, all  satisfying  the  nonpinching  criterion  of  Eq. 

(2).  The  experiments  include  those  of  Sandia  Labora- 
tories, 1 Cornell  University,2  and  Physics  International 
high2  and  low  voltage. 1 Beam  > ranges  from  1.78  to 
4.  91,  v/y  from  0.  9 to  9.  9,  mirror  ratios  up  to  6.7, 

B,  from  2.  5 to  9 kG,  and  the  observed  transmitted 
fraction  from  30%  to  96%,.  We  will  calculate  the  trans- 
mitted fraction  from  Eq.  (7)  if  the  right-hand  side  of 
Eq.  (13)  is  negative.  Otherwise,  the  transmitted  frac- 
tion will  be  the  sum  of  the  core  area  fraction  from  Eq. 
(13)  and  the  adiabatic  contribution  outside  the  core,  ob- 
tainable by  integrating  Eq.  (6)  from  the  core  radius 
to  »V 

Figure  2 shows  a comparison  with  data  from  Ref.  4, 
where  transmitted  fractions  through  20-  and  100-cm 
plasma  channels  are  shown  by  triangles  and  circles, 
respectively,  and  the  calculated  results  are  shown  as 
the  solid  curve.  The  calculation  agrees  well  with  the 
data  of  both  20-  and  100-cm  transport  for  mirror  ratios 
less  than  5.  For  .»/  = 5-6,  the  calculation  deviates  from 
the  100-cm  transport  data.  This  can  be  understood  by 
noting  that  the  total  traveling  time  for  the  reflected 
beam,  which  can  be  roughly  calculated  by  using  the 
number  of  reflections  in  Eq.  (9),  becomes  comparable 
with  the  beam  lifetime  (-100  nsec)  for  such  cases 
(100  cm,  M=5— 6).  Hence,  electrostatic  reflection  in 


the  diode  must  have  been  less  than  100%,  efficient.  Also 
shown  in  Fig.  2 as  a dashed  curve  is  the  adiabatic  trans- 
mission from  Eq.  (7).  Except  for  large  mirror  ratio 
(Af  = 5—6)  and  long  plasma  channel  (100  cm),  the 
adiabatic  transmission  alone  cannot  account  for  the  ob- 
served transmission. 

Table  I shows  a comparison  with  data  from  Refs. 

1-3.  Experimental  beam  and  plasma  parameters  are 
shown  at  the  top  half  of  the  table  and  results  from  the 
calculation  are  at  the  lower  half.  For  Ref.  1,  the  trans- 
mission was  entirely  adiabatic  [Eq.  (7)]  and  agreement 
with  the  data  is  extremely  good.  For  Ref.  3,  one-third 
of  the  calculated  transmission  was  adiabatic  and  two- 
thirds  was  due  to  enhancement;  agreement  with  data  is 
also  good.  For  Ref.  2,  adiabatic  transmission  and  en- 
hancement are  about  equal,  and  the  calculation  falls  at 
the  lower  edge  of  the  data  range.  The  range  (83— 96% ) 
was  the  result  of  various  gas  pressures  (150—900  p) 
in  the  guiding  channel. 

In  summary,  our  model  explains  well  the  results  of 
several  experiments  with  wide  ranges  of  beam  param- 
eters and  transmission  fractions.  We  acknowledge  the 
presence  of  other  potentially  important  mechanisms, 
such  as  formation  of  a virtual  cathode2  in  cases  of  low 
transport  chamber  gas  pressure,  beam-plasma  inter- 
actions which  can  contribute  to  beam  transport  losses, 
beam  diamagnetic  effects,  and  time  variation  of  beam 
parameters.  However,  the  good  agreement  of  our  cal- 
culation with  the  experiments  strongly  suggests  that  the 
mechanism  we  propose  is  often  dominant  in  enhancing 
the  beam  transmission  above  the  adiabatic  prediction. 
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Naval  Surface  Weapons  Center 
ATTN:  Code  WR43 

ATTN:  Code  WA501 


SAMSO/IN 

ATTN:  IND,  Maj  D.  Muskin 

SAMSO/MN 

ATTN:  MNNH 
ATTN:  ffti 

SAMSO/SK 

ATTN:  SKF,  P.  Stadler 

Conmander  in  Chief 
Strategic  Air  Command 

ATTN:  NRI-ST1NF0,  Library 


i 


DE PARTMENT  OF  ENERGY 


Department  of  Energy 
Albuquerque  Operations  Office 

ATTN:  Ooc.  Con.  for  Tech.  Library 

Department  of  Energy,  Library,  Rm.  G-042 

ATTN:  Ooc.  Con.  for  Class.  Tech.  Lib. 

Department  of  Energy 
Nevada  Operations  Office 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

University  of  California 
Lawrence  Livermore  Laboratory 
ATTN:  J.  Nuckolls 
ATTN:  L-153 
ATTN:  L-18 
ATTN:  L.  Woodruff 

ATTN:  Doc.  Con.  for  Tech.  Info.  Dept.  Lib. 

Los  Alamos  Scientific  Laboratory 

ATTN:  Doc.  Con.  for  Reports  Lib. 

Sandia  Laboratories 
Livermore  Laboratory 

ATTN:  Doc.  Con.  for  Library  & Security 
Classification  Division 

Sandia  Laboratories 

ATTN:  Doc.  Con.  for  L.  Vortman 

ATTN:  Doc.  Con.  for  3141 

ATTN:  Doc.  Con.  for  5240,  G.  Yonas 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp. 

2 cy  ATTN:  Tech.  Info.  Services 
ATTN:  P.  Mathur 

Agbabian  Associates 

ATTN:  M.  Agbabian 

Applied  Theory,  Inc. 

2 cy  ATTN:  J.  Trul io 

Avco  Research  A Systems  Group 
ATTN:  Research  Lib. 

Battel le  Memorial  Institute 
ATTN:  R.  Klingsmith 

BDM  Corp. 

ATTN:  Corporate  Library 
ATTN:  A.  Lavagnino 

BDM  Corp. 

ATTN:  R.  Hensley 

Boeing  Co. 

ATTN:  Aerospace  Library 
ATTN:  R.  Carlson 

California  Research  4 Technology,  Inc. 

ATTN:  K.  Kreyenhagen 
ATTN:  Technical  Library 

Dlkewood  Industries,  Inc. 

ATTN:  L.  Davis 

EG4G  Washington  Analytical  Services  Center,  Inc. 
ATTN:  Library 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Ford  Aerospace  4 Comnunications  Corp. 

ATTN:  Tech.  Info.  Services 

Ford  Aerospace  4 Communications  Corp. 

ATTN:  Technical  Library 
ATTN:  MS  G30,  D.  McMorrow 

Gard,  Inc. 

ATTN:  G.  Neidhardt 

General  Electric  Co. 

Space  Division 

ATTN:  J.  Peden 

General  Electric  Company-TEMPO 
Center  for  Advanced  Studies 
ATTN:  DAS  I AC 

I IT  Research  Institute 

ATTN:  Documents  Library 

Institute  for  Defense  Analyses 
ATTN:  IDA,  Librarian 

Ion  Physics  Corp. 

ATTN:  H.  Milde 

IRT  Corp. 

ATTN:  R.  Mertz 

JAYCOR 

ATTN:  E.  Wenaas 

JAYCOR 

ATTN:  R.  Sullivan 
Kaman  AviDyne 
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